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Abstract 

Exploring a classical solution of the non-linear sigma model for a quantum 
Hall ferromagnet, a skyrmion- magnon effective hamiltonian is obtained via 
the collective coordinates method. Using the Feynman- Vernon functional 
integral formalism for this model we find the temperature dependent transport 
coefficients which characterize a single skyrmion dynamics. 

I. INTRODUCTION 

The concept of skyrmion was introduced by T. Skyrme when he represented the nucleon, 
a fermion, as a topological excitation (soliton) of the bosonic pion field jlj. In the quantum 
Hall ferromagnet, skyrmions are the solitons in the spin space of the non linear field theory 
for electrons, and have spin much larger than |. 

The theory for the skyrmions, in the 2-D electron gas subject to strong magnetic 
fields, was developed by Shondi et al || based on an effective generalized non linear sigma 
model (0(3)) which was deduced from the Landau-Ginzburg theory of spinful electrons ||. 
Shondi's conjecture was experimentally confirmed by Barrett and co-workers in the optically 
pumped nuclear magnetic resonance experiments J3J. From both works, || and ||], it can be 
seen that this topological excitations are detectable in the quantum Hall ferromagnet when 
the exchange interaction dominates over the Zeeman coupling. This situation is possible due 
to the smallness of the effective Lande factor (g) in some particular materials, as in the GaAs 
case. If this is so, the simple creation or destruction of an electron in the fully occupied 
spin polarized Landau level, v = 1 for instance, do not produce the low energy configura- 
tion to be considered the lowest lying charged excitation. On the contrary, a non-uniform 
cylindrically symmetric spin configuration is produced, in which the spin at the centre is 
pointing in the wrong direction and smootly rotates as one moves radially outwards until it 
reaches the right direction according to the polarized state of the 2-dimensional electron gas. 
Later on, Hartree-Fock numerical calculations || showed that the energy of the quantum 
Hall skyrmion excitation is always smaller than the excitation energy of the localized spin | 



*E-mail: caldeira@ifi.unicamp.br 



1 



quasiparticle and quasihole, supporting the idea that the skyrmions are for sure the relevant 
quasiparticles in the g — > limit. 

In practice it is convenient to regard the skyrmion as a topological object. From this 
point of view they can be described as the non-singular mapping of the compactified 2-D real 
space (Sf^) onto the 2-D internal space (Sl nt ) whose elements are the unitarily normalized 
spin vectors which cover a 2-D spherical surface. All non-singular maps of 5*2 onto S2 can 
be classified into homotopy classes that form a group isomorphic to the group of integers. 
And, it is using this set of integers, that we define the topological charge or Pontryaguin 
index Q, that classifies the skyrmions. This topological charge can be shown to be a simple 
measure of the degree of mapping of one surface onto the other and, as a conserved number, 
helps in finding explicit expressions for the spin configurations (see H, for instance). 

Until now many theoretical studies were devoted to the properties of the quantum Hall 
skyrmions. Most of them are done following two main approches; the quantum compu- 
tational calculations and the non-linear semi-classical field theory. In the first case, the 
quantum Hartree-Fock calculations are shown to be in excellent agreement with exact diag- 
onalizations for small skyrmion sizes, as is demostrated in H. However, when hydrostatic 
pressure is applied to the quantum Hall ferromagnet the g factor almost vanishes. In that 
case there is strong experimental evidence |J that skyrmions of large size (about 33 reversed 
spins) are dominant. In this particular situation computational calculations are not appro- 
priate, making the development of new approaches in the semi-classical field theoretical 
methods the natural way out to explain the experimental data. 

In this paper, starting from a specific Lagragian density 0, we investigate the influence 
of the ferromagnetic spin wave modes in the temperature dependent dynamical properties 
of an isolated skyrmion. Our model basically consists of an effective field theory in which 
the single skyrmion momentum is coupled to the magnon bath total momentum. To begin 
with, sec II is devoted to obtaining the ferromagnetic spin wave modes by computing the 
quantum fluctuations of the quantum Hall ferromagnet around the non-linear sigma model 
classical solutions. Later on, using the standard collective coordinate method a second quan- 
tized hamiltonian that simulates the scattering of the spin waves by the skyrmion in 2-D is 
obtained. 

In sec. Ill, using the well known Feynman- Vernon formalism we obtain an effective 
equation of motion for the centre of mass of the skyrmion. This equation, that contains all 
the information about the skyrmion-magnon interaction, is written in terms of a temperature 
dependent damping matrix function. In our treatment the possible contribution of inelastic 
process (Cherenkov like radiation) to the soliton dynamics is considered. Finally sees. IV 
and V are devoted to the explicit calculation of the damping constant matrix elements that 
define the skyrmion mobility and to present our conclusions. 
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II. EFFECTIVE SKYRMION-MAGNON HAMILTONIAN 



As we said, the quantum Hall skyrmions can be described in terms of a unit vector field 
associated to the local pseudo-spin orientation as was proposed by Shondi [0. The effective 
lagrangian density of the system can be written as 



where 



£(n) = aA [n(r)] • d t n{r) - V[n] 



A [n(r)] = (d t n x d T n) 



(1) 



(2) 



is the vector potential of the unit monopole, and the effective potential functional is that of 
a ferromagnet with long range interaction arising from the Coulomb repulsion V(r) = e 2 /r 
between the underlying electrons. Explicitly, 



1 r 

V[n] = a [Vn(r)] 2 + gpp B n{r) • B + - / d 2 r'V{r - r')q(r)q(r') 



(3) 



where the deviation from the uniform background density p = 1/2tcIq is the skyrmion density 
q(r), which in terms of the n field has the following form 



q(r) 



57T 



£„j,n ■ (dfjii x d u rv) 



(4) 



Here a' is the spin stiffness and lo stands for the magnetic lenght. 

To take advantage of the internal space spherical symmetry, it is convenient to define 
the n vector field by means of two escalar fields G(r, 9) and $(r, 9) such that, 



n 



n" 



n 



sin[e(r, 9)] cos[$(r, 9)) 
Bin[e(r, 9)] mn[$(r, 9)] 
cos[6(r,0)]. 



(5) 



In terms of these new fields the effective potential functional (|3|) normalized to the Zeeman 
energy reads, 



V = a' 



(ve) 2 - sm 2 e (v$) 2 l + gnBp{i - cos e) 



e 

+—q(r) 



d(e, $) 



d(r, 9) 



q{r 



d(0, $) 



d{r',9') 



dr 



r> 



(6) 



where 



9(e,$) 

d(r,9) 



stands for the Jacobian of the transformation and the skyrmion density is 



simply q(r) = sin[0]/47rr. 



From (||), the localized static solutions in spin space (skyrmions) can be derived from 
the variational equations 5V/5Q = = 0. Explicitly, 



= a' 



'n^ sin 26 /T -^ x2 
V 2 (V$) 



gpsBpsin 6 + e 2 q(r) / q(r' 



<9($, |r 



„/r— l-s 



d(r, 9) 



dr 



(7) 
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= a' sin 2 0V 2 $ + sin20(V$).(V0) + e 2 q(r) q{r') 



d{\v 



d{r, 9) 



dr' 2 . 



(8) 



This set of equations with the boundary condition O — > as r — > oo and cylindrical sym- 
metry was already solved by numerical methods ||. However, we will use a variational 
approach, as in , in which, starting from the non-linear sigma model (NLcr) solutions [|J 



= arccos 



~2Q _ 4 ^2Q ^| 
-2Q _|_ 4^2Q J 



and 



$ = arctan{(-) Q }, 



(9) 



(10) 



an optimum value of A is determined by balancing the Coulomb and the Zeeman contribu- 
tions in (Bf). For the topological charge Q = J q(r)d 2 r = 1, A has the value 



A = 0.558/ (£|ln£|H 
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gfi B B 



(11) 



where g is the reduced Lande factor which varies from 0.02 in the normal GaAs material to 
0.002 when an hydrostatic pressure is applaied ||. 

Therefore, at T = and near v — 1, we will assume that the quantum Hall ferromagnet 
is a non- interacting skyrmion gas given by (If)- (pp. However, when T^0 the ferromagnetic 
spin wave modes of the system are excited. This kind of quantum fluctuations are given 
essentially by the functional expansion of the effective potential V[0, $] around the NLa 
model solution (|9|) and ([TO) up to second order. As the optimized size for the skyrmion gives 
a negligible first functional derivative contribution, the expansion generates the following 
equation for the fluctuations of the field, 



{a'V 2 + U[G, $]}w m , k (r) = E m (k)u mM (r), 



(12) 



where 



with 



U[e,$\ = a' cos(2<d)(V$) 2 + gfi B Bp cos Q + e 2 q(r) J q(r')K[Q,^}dr' 2 (13) 



K 



<9 $(r), 



9([r-r , r 1 ,9(rQ) 



d(r',e') 



d{r, 9) 



(14) 



Substituting @ and ( [10]) in ( |T3| ) the following Schrodinger-like equation, for Q — 1, is 
obtained 



{-V 2 + C/(r)}« m>k (r) = e m M mik (r) 



(15) 



4 



where 



and 



E m - 9 ^ (16) 



«7| 2VA% 2 2gfi B BX 2 
{T> r 2 (r 2 + 4A 2 ) 2 7r(r 2 + 4A 2 )' 1 ' 

Notice that due to the rotational invariance of (§) and the radial symmetry of (fL0|), the 
Coulomb term does not contribute at all to the fluctuations within this approximation, so, 
far from the scattering centre, the fluctuations will be mainly dominated by the potential 
generated by the skyrmion itself and by the presence of the magnetic field. 

Despite of the complex form of (|i~7|), the solutions of ([L5|) for the m-th cylindrical wave 
component of the ferromagnetic spin waves confined in a cylindrical region of radius L can 



be written as the sum of two contributions [10] 

Umn oc \{H^{k mn ry m * + J2e- 2t5 ^H^(k ln r)e^}, (18) 
1 i 

where H^' 2 \kr) are the Hankel functions of first and second kind respectively and 5 m i is 
the phase shift matrix which connects the m and I angular momentum channels. The first 
term of the right hand side of ([18]) correspond to the incident cylindrically symmetric wave 
whereas the second term is associated to the shifted Z-th component after the interaction 
with the potential U(r) has taken place. As can be seen from (fl7|) , the potential involved in 
the scattering processes is cylindrically symmetric, therefore 5 m j is a diagonal matrix and 
(|i~8l) becomes 

u mn cx \{H\2\{k mn r) + e- 2lS ™H^(k mn r)}e im t (19) 
To calculate the phase shifts for each angular momentum component we used the Fred- 



holm method ||11|| . In so doing it is convenient to start from the well-known relation, 

irA{E) cot 5{E) = l + vJ^^dE', (20) 

o 

where V means the Cauchy principal value, and A(E) can be obtained in first approximation 

as 

A(E) = -(E\U\E). (21) 

In computing A(E) we will use a cylindrical basis solution of the non perturbed Schrodinger 
problem with circular boundary of radius L. When L is large enough, the 2-D cylindrical 
basis can be written as 

/2/\/f \ 3 

HE>=(^r) J m (kr)e- e . (22) 
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However looking at (|I7D it can be seen that, as r — » 0, U(r) — > oo and hence every 
component of fl2"2"D must vanish at r = 0. This condition is automatically satisfied by the 
m^O J m (kr) functions. Therefore substituting ( p2|) into fl2T| ) we get 



Am{x) = JL \l m {x)[K m+x {x) + K m _x{xj\ - K m (x)[I m+ i(x) + I m _i(x)] - (4mx)~ j 



with 



8V2vr 

+2e z I m (x)K m (x) } (23) 



e 2 /^o „ _ QHbB ^ 



where we have used the value of the spin stiffness e 2 /32^/2^ll of the quantum Hall ferromag- 
net. On the other hand I m and K m are the modified Bessel fuctions and x = 2k\. Then, 
the phase shift for each m angular momentum component is given by 

6m = arctan (iT^SS)) ' (25) 

where we have used A(2Xk) to denote the expression 

OO A ( \ 







So, until now, our model for the quantum Hall ferromagnet is composed by a set of 
static cylindrically symmetric skyrmion excitations given by (fl)-(|l]]) and non interacting 
spin waves modes described by ([19]) and (F23|)-(|26D with frequencies 



Umn ~^M + "ft - ' (27) 

where the momenta k mn can be determined from the boundary condition u mn (kL) = 0. 
Therefore using the asymptotic expansion of (|T9|) we get 



fc m „ = ^ + — (28) 
2L 2L L K ' 

On the classical level, due to the translational invariance of the system, the skyrmions 
can move freely, and the position of its centre of mass will be enough to describe its motion. 
However, if we look at the problem from the quantum field theoretical point of view, the 
skyrmion will be a particle and its centre of mass a dynamical variable. In this picture, due 
to the excitation of the spin wave modes at finite T, not all of the degrees of freedom of the 
system will contribute to the skyrmion formation, this can be interpreted as the origin of the 
residual skyrmion-magnon interaction that leads to a non trivial dynamics of these objects. 
To show that, we explicitly evaluate the equation of motion for the B field (our skyrmion 
field now) from the equations ([[]), ([§) and (^). Considering cylindrical symmetry for non 
static skyrmions and neglecting higher order spatial derivatives, the effective equation of 
motion for 0(r,t) can be written as, 



6 



4 e +2a' 



v 2 e 



sin 26 



gpHsB sin = 0, 



(29) 



where c = a' /a. Comparing (|7|) and ( |29]) it can be seen that to describe the skyrmion 
dynamics interacting with a magnon thermal bath in a simple way, we can start from the 
one field Lagrangian density associated with the z pseudo-spin component which reads 



where 



v[e] 



a 



C 



(V0) 5 



sin 2 9 



+ gfi B Bp(l — cos 6). 



(30) 



(31) 



Therefore, our assumptions have simplified the initial problem given by ([!])- ([|) into an 
effective two dimensional one escalar field theory that describes the dynamical skyrmion and 
the possible scattering states given by ferromagnetic spin waves. 



At this point, our model to treat the transport properties of the skyrmions is given by 
(j30l) and (|31~1) . The correct quantization of ( |30D can be done using the canonical formalism 
via the method of collective coordinates [|l^]. Doing this we approprietely treat the zero 
frequency mode associated to the translational invariance, and the second quantized version 
of H (obtained from fl30D and ([H])) which describes the momentum-momentum coupling 
between skyrmions and magnons (see |13[ for details) can be written as 



H 



2M 



mg 



~)~ ^ ^ ^^mnb mn b r< 



(32) 



In (|32"D, P s and P m3 denotes the skyrmion and magnon momentum operators and u v 
is given by Q2"7|). Explicitly, 



■mg 



^mn,klbmJ ) kl + ^ ^mn,kl (bmnbkl ~ b mn b 



t 

kl) > 



(33) 



mn,kl 



mn,kl 



with 



D 



mn,kl 



LOr, 



(Jj„ 



G 



mn,kl t 



(34) 



^m.n 



OJr. 



G 



rnn,kU 



where G m n,ki is a 2-D vector given by 

GmnM = J u kl (r)Vu mn {r)d*r 

and 



(35) 



(36) 
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M = J<yu (r)) 2 d 2 r, 



(37) 



where u denote the localized solution (§) with Q = 1 whose gradient is associated to the 
zero frequency mode. 

The relations (|32"D-(|37D are a generalization for the 2-D case of the results obtained in 
13fl . Notice from ([33|), that the magnon momentum operator consists of two contributions. 



A block diagonal term, which only couples excitations with different fc, as can be seen from 
the antisymmetric behaviour of D mn ^ under inversion of index, and that is responsible for 
the low energy scattering of magnons by the skyrmions. This term commutes with the total 
number of magnons, so any kind of transition involving processes which do not conserve the 
number of magnons is not allowed. Actually, this contribution remind us of the problem of 
scattering by a hard (but not fixed) object because there are no internal degrees of freedom 
in our model. On the other hand, the block off-diagonal terms in ( |3"3"D describes the emission 
or absorption of magnons by the skyrmion. In some specific situations, during the Hall 
measurements or in the high temperature limit, for instance, the effective motion of the 
skyrmion can be affected by the off-diagonal contribution. Although we do not attempt 
to study the contribution of the inelastic process to the skyrmion dynamics in this paper, 
we will show that in the long time approximation this term do not contribute at all to the 
transport properties. 



III. THE FEYNMAN- VERNON FORMALISM 



Based on the functional integral method of the Feynman- Vernon formalism fT| , Castro 
Neto and Caldeira developed a new model for dissipation in quantum mechanics [[0| to 
study the polaron dynamics [13| without appealling to the kinetic theory. In this work we 
will follow the same kind of approach, but now in the 2-D physical space, in our study of 
the effective skyrmion dynamics in the interacting system (skyrmion plus magnons). 

The key point of the method is always to compute the reduced density operator (p s ) for 
the particle of interest (the skyrmion) by tracing the mesons (magnons) coordinates out of 
the problem. The density operator of the interacting system evolves in time according to 

pit) = e-^p(0)e^, (38) 
where H is given by (|32]) and p(0) is assumed to be decoupled at t = as 

p(0) = p s (0)p m ,(0), (39) 
where s and mg denotes skyrmion and magnon respectively and 

XV (c~^^ m9 ) 

As it was shown in ||13|| , the reduced density operator in the coordinate representation 
for the particle of interest (skyrmion) can be written as 



S 



p(x, y,t) = J J dx'rfy'j(x, y, t; x', y', 0)p.(x', y', 0), (41) 

where J is the skyrmion superpropagator which can be explicitly written in terms of the 
influence functional .F[x, y] as 

x y 

J = Jv*J Pyes( s °M- s °M):F[x,y], (42) 

x' y> 

where the action associated to the free skyrmion motion is 

S [x] = J M^dt, (43) 
o 

with M given by (]37f) , and the influence funtional in the coherent state representation is 
written as 



a a 



& 7* 
e (a*(0)a+a*a(t)+ 7 (0)a*+a 7 *(t)) j 

Here we have introduced Si as the interacting skyrmion-magnon action, and p mg stands for 
the magnon density operator in the coherent state representation, 



/iTi 
{ y J2 a *mn V -On Q* mn - £ flU., 
II ran ran 

+HM X • ^ ^>mn,kl a *mn a kl 



mn^ran^ran 



rnn,kl 



+ hM X • Cmn,kl( a rrm®kl + ®*mn a kl) } dt '> ( 45 ) 



ran,kl 



and 



o* — - — BHuJmr, \l3mn\ 2 \~fmn\ 2 

^,7) =n — (tt^h — • ( 46 ) 

ran v x ° / 



Notice from (fig) that changing to the Lagrangian formalism through the coherent state 
representation we have simplified the problem (|32"D, obtaining an interacting skyrmion- 
magnon action quadratic in a. Therefore the equation of motion obtained from the varia- 
tional problems SSi/Sa = 5Si/5a* = can be formally solved. 

The linear Euler-Lagrange equations obtained from (051) can be written as 



iM ±-Y^ ~Dmn,kl a kl + 0, (47) 

ki 
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a* mn -iu mn a* mn + iMx^ T> mnM a* kl + 2C mnM a u = 0, (48) 

kl 

whose solutions, in terms of the functional operator W nm are of the form, 



(r)a r 



<{r) = e* 



<+ E (Wmn,w(T)a fcI )* e - <w « t 



(49) 



(50) 



The functional operator introduced above obeys the integral equation 



WWx(r)] = / W°(r')dr 



with 



+ / E ^ mnm {r')W pqM {r')dr' (51) 

PQ^ran 



W° mnM = iM x ■ {D mr ^e- 4 (— + 2C* mnM e^- + ^ t } . (52) 



Now, using (H|) and (0) we are able to solve the functional integrals in (|4]) in the 
stationary phase approximation. After evaluating some integrals, the influence functional 
can be further expressed as 

^■yi= t/l 4 r — v (53) 

Qet ^1 — 1\ mn L mn.kl J 

where T mnM is given by 

r mn ,ki = W mnM [y] + Wl nM [x\ + W mnm [y]Wl qM [x\. (54) 
and the ocupation number for magnons is, 

Nmn exp((3?ko mn ) - 1 ' ^ 

To go on in computing the influence functional we need to explicitly solve fl5"l"|). This can 
be done up to any order by iteration, however as this relation is nothing but the amplitude of 
scattering from one state mn to another kl through virtual transitions involving the state pq, 
we hope that in describing the low velocity skyrmion dynamics only small energy processes 
are involved and therefore a second order iteration will be sufficient to describe the main 
features of the transport properties of the system. Doing this, and after some algebra, the 
influence functional can be expressed as 
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F[x,y] = e^ I +^e-^ R +* R \ 
Where the following definitions were used, 



(56) 



t t 



= jdtjdfY, {[xdt')+y t (t')}rK(t'-t")[x 3 (t")- yj (t")}}, (57) 







t t 



$ / = J dt'f dt" j2 {h (*')+ ^ (*')] F L(*' - o h (O 

o o 'J- 1 



v,- (*")] } , (58) 



= / ^ j dt" £ - O Vj (O- ^ (^)A£(t' - O 2/, (O] (59) 



o o hi= l 



t t 



¥ = J dt' J dt" J2 [*i (t")^lj(t' - 1") v 3 (0+ x t (t')A^(t' - 1") vj (f) 

o »J' =1 



(60) 



where and Ajj are 2x2 matrices with elements involving the components of the 2-D 



vectors ~D m n,ki and C m „ )fc j g iven by (|3^)-(Pq). Explicitly 



(^mn + + 2N mn N kl )(u mn + l^y) 1 



mn,kl 



2u) mn U) k \ 



MiJmn, kl) cos(uj k i - u mn )t, (61) 



r i,iW = 2^ — ~ Mij{mn, kl) sm(w u - u mn )t, 



mn,kl 



2Ld rnn UJkl 



(62) 



a Ft V™^ V ly inn 

A m'W = E ^ — ~ Mi,j(mn,kl)cos(uj k i + uj mn )t, (63) 



mn,kl 



2oj mn u) k i 



\ I U\ (^mn + N k l + 2N mn N k i)(uJ mn — LOklf , . . 

A i iW = E o -Mj(mn,A;Z)cos(w w + u; mn )t, (64) 



mn,kl 



2u mn u k i 



i r (i) 

I ran 

l^mn,fcZ 



•^>i = n ,^(2) l2 (65) 



where G^^ kl are the components of the G vector given by (|36|) . Finally, substituting (|56| 
in (H2I) the skyrmion superpropagator becomes 
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x y _ 

S[x.,y]+ih(0 R -^ R ) 



J= The -py e *[**\-W-*~). (66) 



r 



Now, after having traced the magnons coordinates out, the reduced skyrmion action 
contains all the information of the interaction with the thermal bath and can be written in 
a simple way as 

S[x,y] = 5 [x] - S [y] + ih{& + ¥). (67) 



Before analyzing effective skyrmion motion it is convenient to rewrite fl67D in terms of a 
new set of variables, R = (x + y)/2 and r = (x — y). The equations of motion associated 
to the action 5*[R, r] can be written as 



t 



Ri (t) + J ji,j(t - Rj (t')dt' = 0, (68) 



o 



r, t (t) - J j i:j (t - t') rj (t')dt' = 0, (69) 
o 

where we have introduced the matrix damping function, 7»j(t — t'), that characterize in a 
general form the mobility of the 2-D solitons. In terms of fl62"D , (0) and (|65|) the matrix jij 
is given by 

7«j(< - = ~~ - + A^(i - t')} . (70) 



The equations of motion (|68| ) and ( ]69| ) are generalizations of those obtained in the case of 
the soliton dynamics |15| describing a quantum Brownian motion with memory effect. The 
main differences are that the present treatment allows for the study of the 2-D topological 
excitations and that the Cherenkov term accounts for the possibility of mesons (magnons) 
emission or absorption by the soliton (skyrmion). This contribution can be important, even 
for small velocities, if we are studying the time evolution of the soliton. 



On the other hand, as in [0, the real part of the exponent in (]67|) is related to the 
matrix diffusion function 

A.(t) = ^(r£ + A5), (7i) 

that can be computed in close analogy to the matrix damping functions as shown in the 
next section. 

To simplify the analysis of ( [70]) the explicit expression for the damping function will be 
written as a sum of two terms, 

7w(*-0 = 7S(*-0 + 7&(*-0, (72) 

where 
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- N kl )(uJ mn + UJ kl ) 2 (uj H - UJ 

ran I , , / , 7 \ / \ / . . ;\ 

Hi = 77 L o Mij{mn, kl) cos(uo kl - u mn )(t -t), 

ranM 2uJ mn UJ kl 



(73) 



and 



! _ 2h _ (N mn + N kl + 2N mn N kl ) 

(u mn - to k i) 2 (uj k i + u mn ) cos(u k i + UJ mn ){t - t'). (74) 

As it can be seen from (|36|) and (|65|) the matrix elemements A^ij essentially involve the 
single particle momentum operator between two different ferromagnetic spin wave states 



given by (|19D . Therefore in close analogy to the spectral fuctions introduced in |16| we will 
define the 2x2 matrix scattering function, Sij, as 

Sij{u, u') = M itj (mn,kl)S(uj -u mn )S(uj' - u M ), (75) 

mn,kl 

that allow us to rewrite the expressions for 7/ • and jfj in the following form 

n 



iF/t) = -0(t) J duj dcSSijfautyNiu,) - N{lS)) 

^— — ^ >- cos(cj' - u)t, (76) 

which corresponds to the elastic scattering process with purely instantaneous memory when 
we assume long time approximation, as in ||13||. On the other hand 



UJUJ 

\N{uj) + N{uo') + 2N{uj)N{uj')} cos(cj + uj% (77) 

is related to the Cherenkov radiation and only contributes to the soliton dynamics when 
transitions between states with different energies are taken into account. So, we will refer 
to this contribution as inelastic. In Q7T3| ) and ( [77] ) 0(t) ensures the causality principle. 



IV. THE DAMPING MATRIX 

In the first place let us analyze the J* At) contribution to the damped motion of the 
skyrmion. Notice from (|77|) that for long times the cos(o/ + u)t term oscillates rapidly giv- 
ing no net contribution to the damping function matrix elements. This is true unless uo' + uj 
is close to zero, but at the same time the factor (uj — uj') 2 prevents this situation. So, up to 
first order in the solution of ([51]), the inelastic contribution 7/ ? -(i) to the damping function 
can be neglected. However, if we consider processes with characteristic times of the same 
order as (uo' + uo)~ l , when the long time approximation is no longer valid, this term will play 
an important role in the soliton dynamics. 
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Now, the analysis of (|76|) becomes simpler if we change the frequency variables to ip = 
(u/ + u)/2 and r\ = (a/ — uj). In this way 



oo oo 



75(0 = %W)j*p f d vSi M vWf(<P, ,) NiV + " /2) " N(V - " /2) (78) 

-oo ^ 

where 

/(*,,) = 2 + ^§ + <^§. (79) 
^ - ?7/2 9? + 77/2 

As in ([T7]), the expression f[78|) vanishes for long times due to the rapid oscillations of 
the cos(yyt) term unless r\ is close to zero. If this is so f(<p,r]) can be taken as constant and 
we can assume that 

N(ip + r,/2)-N(ip-r,/2) ^ dN(<p) 

— , (m 

r] dip 

therefore ( [78f) can be rewritten as 

oo oo -— 

iUt) = ^e(t) fdtp j dvAM^Mcosivt), (si) 

-oo ^ 

where we have introduced Ai j(cp) given by 

A id {<p) = limS i>j {ip, V )r ] 2 . (82) 



After integrating fl8~l|) in dr) and using fl55|) we get 

7 £(t)=7 M (T)5(t), (83) 

where 7 y i j(T) is the temperature dependent matrix damping parameter that characterizes 
the soliton mobility in the 2-D case and 8(t) is the Dirac delta function. As usual the 
we will refer to the form (|83|) as the Markovian approximation, because this results have 
only instantaneus memory not depending on the previous motion of the particle. Explicitly 

%j( T ) is s iven b y 

%A T ) = mJ (1 _ ePi^B+n^ A^. (84) 
o v ' 

In order to compute the temperature dependence of % ;(T) we need to explicitly evaluate 
the scattering function Sij(<p, if). The first step in doing so is to get an analitic expression for 
Sij(u),u'), and therefore the matrix elements of ( pop must be evaluated. As it can be seen 
from (|l9l) the x and y components of the G vector defined by (|65|) are equal, in agreement 
with the isotropic character of our model. Therefore the function matrix A4 m n,M can be 
written as M. = G mn ^{\- where 1 is the unitary matrix. Then we need to compute only one 
of its componenets, for instance, 
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f J d sine d 

G mnM = I rdr I d9u* kl (cos9- — )v, 



dr r 89' 



(85) 



o o 

which, after some integration and assuming large values of L becomes 



G 



mklk' 



A Wp k _L ^UkV)K1, _ U\ _ 4„U\WxfU , fc/M _ a( 4 )~ ^ 



,( 4 ) r^t , ^a( 6 )" 



where C;(/c) stands for 



Ci(fc) 



cos fct 

t 



and the set of expressions 

A (n) 

are defined by 

A (1) = ( i+m ) e 2*<5m _|_ e 2i«5, e -if (Z+m) 



if Ag + - fc')A^ - C,(* + fcOAg (86) 



(87) 



(88) 



A (2) = e if (i-m) _ e -q(l-m) e 2i(5 t -5 m ) 



(89) 



A (3) = e if (l+m) e 2iS m + e 2i<5 ie -i|(/+m)^ 



(90) 



A (4) = e tf (I-m) + e -tf (l-m)g«(«,-*») 



(91) 



A (5) = A (2) + e 2i5 ig -if (i+m) _ e -2i<5 me if (J+m) 



(92) 



A (6) _ gif («+m) e -2i5 m + e 2i<5 ie -if (I+m) 



(93) 



Expression (^3) allows us to write a continuum version of Q73 ) for the matrix damping 
function diagonal elements S^. Explicitly 



L 4 /" r 



(94) 



where 



2M' 



(95) 



Now, using (p6|)-(pOD we can evaluate (Q) and obtain an explicit expression of Si t i(u,u}') 
which, by the time, allow us to get the analytic form of A iti in terms of the variables <p and 
r) introduced before. Then, the limit given by yields 



15 



AM = (^jf) 2 G&)<P 2 (96) 

with 

oo 

£(y?) = 2£sin 2 (5 m+1 -5 m )- (97) 

m=l 

where the phase shifts 5 m are also functions of ip. 

The matrix damping coefficient elements ( pip can be written now as 

~ . = (*™\ e ^B 7 gv^gw , (98) 



In computing this final expression (|9q) we need a numerical evaluation of the phase 
shift functions defined by (p3"D-(p6|). This was done taking the magnetic field B = 9T for 
three different values of the Lande factor g or skyrmion size. However, the analysis of the 
high temperature regime of the skyrmion mobility can be done from ( P5| ) without explicit 
evaluation of Q, giving 



that increases linearly with temperature independently of the explicit form of Q(<p)- 

The behaviour of 7jj(T) in units of 70 = ir 2 h 2 /32M 2 A 4 for any temperature, is shown in 
Fig.l. for three different values of the Lande factor. Due to the smallness of the Zeeman 
energy (0.13 meV, 0.013 meV and 0.0013 meV) in all cases the scattering states are excited 
leading to a damped motion of the skyrmion. Although we do not have an analytic expression 
for 7jj(T) the numerical result shows that near T = the mobility presents a power law 
behaviour T a , as can be seen from the log-log plot in Fig.2, where a is a strong dependent 
Lande factor coefficient. 



V. CONCLUSIONS 

Here we have succeded in applying the collective coordinate method of quantum field 
theory to the description of the mobility of the skyrmion excitation in a quantum Hall 
ferromagnet as a function of the temperature. We have been able to treat both elastic and 
inelastic effects of the skyrmion-magnon interaction and show that only the former plays a 
major role in the long time approximation for this problem. We have also shown that in 
this limit the mobility of the skyrmion is linear for high temperatures and behaves as T a 
for low temperatures. 

Although we have directly addressed the problem of the quantum Hall ferromagnet, our 
description can be straightforwardly extended to any system which can be mapped into a 
NLcr model or any extension thereof. An investigation of the experimental evidences of our 
findings is now in process. 

Finally AVF wishes to thank Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP) for financial support, whereas AOC kindly acknowledges partial support from 
Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq). 
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FIGURES 



x 10 
x 10"' 
x 10 s 



2.0 



0.0 




kT (meV) 

FIG. 1. Damping coefficient as a function of temperature for different effective Lande factors. 
The solid line correspond to the case in which g = 0.25, the dotted line to g = 0.025, and the 
dashed line to g = 0.0025 
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FIG. 2. The circles correspond to the case in which g = 0.25, the squares to g = 0.025, and 
the triangles to g = 0.0025 respectevely. 
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